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Abstract 

Convective heat transfer can be enhanced by changing flow geometry and/or by enhancing thermal conductivity 
of the fluid. This study proposes simultaneous passive heat transfer enhancement by combining the geometry 
effect utilizing nanofluids inflow in coils. The two nanofluid suspensions examined in this study are: water-AI 2 0 3 
and water-CuO. The flow behavior and heat transfer performance of these nanofluid suspensions in various 
configurations of coiled square tubes, e.g., conical spiral, in-plane spiral, and helical spiral, are investigated and 
compared with those for water flowing in a straight tube. Laminar flow of a Newtonian nanofluid in coils made of 
square cross section tubes is simulated using computational fluid dynamics (CFD)approach, where the nanofluid 
properties are treated as functions of particle volumetric concentration and temperature. The results indicate that 
addition of small amounts of nanoparticles up to 1% improves significantly the heat transfer performance; 
however, further addition tends to deteriorate heat transfer performance. 



Introduction 

Convective heat transfer can be enhanced by active as 
well as passive methods. While the former usually pro- 
vide better enhancement, it requires additional external 
forces and/or equipment which can increase the com- 
plexity, capital, and operating costs of the system. In 
contrast, passive heat transfer enhancement can be 
achieved by changing flow geometry or modifying 
thermo-physical properties of working fluid. Hence, it is 
generally a more desirable approach when compared to 
an active method. In our previous study [1-3] (Sasmito 
AP, Kurnia JC, Mujumdar AS: Numerical evaluation of 
transport phenomena in a T-junction micro-reactor 
with coils of square cross section tubes, submitted), we 
have shown that coiled tubes provide better heat trans- 
fer performance relative to straight tubes under certain 
conditions. In this study, the potential application of 
coiled tubes using nanofluids to improve heat transfer 
performance is investigated. 

Coiled tubes have been known as one of the passive 
heat transfer enhancement techniques in heat and mass 
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transfer applications due to the presence of secondary 
flows which improve heat and mass transfer rates. They 
have been widely used in process industries, e.g., heat 
exchangers and chemical reactors, due to their compact 
design, high heat transfer rate, and ease of manufacture. 
Aside from their industrial applications, studies of the 
transport phenomena in coiled duct have also attracted 
many attention from engineering researchers. The pre- 
sence of secondary flows induced by coil curvature and 
the complex temperature profiles caused by curvature- 
induced torsion are among significant phenomena 
which can be observed in coiled tubes. Numerous 
experimental [4-8] and numerical [1-3,9-13] investiga- 
tions on heat transfer and flow characteristics inside 
coiled tubes have already been reported. Furthermore, 
reviews on the flow and heat transfer characteristics and 
potential application of coiled tubes in process indus- 
tries and heat transfer application can be found in 
[14,15]. 

It is well known that conventional heat transfer fluids 
including water, oil, and ethylene glycol mixtures have 
poor heat transfer rate due to their low thermal conduc- 
tivity. Therefore, over the past decade, extensive 
research have been conducted to improve thermal con- 
ductivity of these fluids by suspending nanoparticles of 
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diverse materials in heat transfer fluids, called nanofluids 
[16]. Modern technology provides opportunities to pro- 
cess and produce particles below 50 nm. It is also 
expected that nanofluids should provide not only higher 
heat transfer rate, but also good stability of the suspen- 
sion by eliminating possible agglomeration and sedimen- 
tation to permit long-term application [17]. To date, 
several experimental (see for example [18-23]) and 
numerical (see for example [24-28]) investigations to 
characterize heat transfer performance of nanofluids 
have been already reported. Choi et al. [18] showed that 
addition of small amounts of less than 1% nanoparticles 
can double the thermal conductivity of working fluids. 
Vajjha et al. [24] showed that heat transfer rate increases 
up to 94% by adding 10% A1 2 0 3 nanofluid and increase 
up to around 89% by adding 6% CuO nanofluid. In 
addition, the comprehensive reference on nanofluids can 
be found in the book of Das et al. [29], while several 
reviews of nanofluids are available in the literature 
[30-42]. 

It has been shown that coiled tubes geometry and 
nanofluids can passively enhanced heat transfer perfor- 
mance. Now, to maximize the advantages of the heat 
transfer enhancement, we propose to combine both 
techniques simultaneously; i.e., employing the combina- 
tion of coiled tubes filled with nanofluids. Therefore, the 
aim of the study presented here is threefold: (i) to inves- 
tigate the heat transfer performance of various config- 
urations of coils of square tubes, e.g., conical spiral, in- 
plane spiral, and helical spiral, relative to the straight 
pipe; (ii) to evaluate simultaneous passive heat transfer 
enhancement-channel geometry and fluid thermo-physi- 
cal properties-in coiled tubes filled with nanofluids; (iii) 
to study the heat performance of two different nano- 
fluids, water- A1 2 0 3 and water-CuO, in coiled tubes at 
various nanoparticle concentrations. The most signifi- 
cant aspect of this study is to determine the potential 
advantages and limitations of heat transfer enhancement 
of coiled of square tubes filled with nanofluids and pro- 
vide design guidelines for their applications through 
mathematical modeling. 

The layout of the article is as follows. First, the mathe- 
matical model is introduced; it comprises conservation 
equations for mass, momentum, and energy. The nano- 
fluid ther mo -physical properties are treated as functions 
of particle volumetric concentration and temperature. 
The mathematical model is then solved numerically uti- 
lizing finite-volume-based CFD software Fluent 6.3.26, 
the User-Defined Function written in C language is used 
extensively to capture the nanofluid properties. The 
model is further validated against experimental data by 
Anoop et al. [19] in terms of heat transfer performance 
for both base-fluid and nanofluid. Fluid flow and heat 
transfer performance of various coiled tube designs filled 



with nanofluids is evaluated in terms of a figure of Merit 
Defined later. Parametric studies for particle concentra- 
tion and nanofluid type are then carried out. Finally, 
conclusions are drawn and possible extensions of the 
study are highlighted. 

Mathematical model 

The physical model (see Figure 1) comprises four tube 
designs, e.g., straight pipe, conical spiral, in-plane spiral, 
and helical spiral, filled with two different nanofluids 
(water- A1 2 0 3 and water-CuO). We assume that the low 
particle volumetric concentration of nanoparticles (less 
than 3%) in the base-fluid makes it behave like a single- 
phase fluid and there is no agglomeration or sedimenta- 
tion which occurs inside the tubes. A constant wall tem- 
perature is prescribed along all sides of the channel wall; 
the nanofluid is assumed incompressible and Newto- 
nian. Furthermore, to ensure fidelity of the comparison 
of heat transfer performance for each tube design, the 
total length of each tube design is kept constant. Since 
this study relates only to laminar flow, a precise numeri- 
cal solution is adequate to simulate reality very closely. 

Governing equations 

In the tube, fluid flow and convective heat transfer are 
taken into consideration. The con-servation equations of 
mass, momentum, and energy are given by [24] 

V • (pnfU) = 0, (1) 

V • (p nf u ® u) = - Vp + V • [/x nf (Vu + (Vu) T )] , (2) 

V • (p nf Cp, nf uT) = V • (fc n fVT). (3) 

In the above equations, p n f is the nanofluid fluid den- 
sity, u is the fluid velocity, p is the pressure, ^ n f is the 
dynamic viscosity of the nanofluid, c p>nf is the specific 
heat of the nanofluid and k n f is thermal conductivity of 
the nanofluid. 

Constitutive relations 

Thermo-physical properties of nanofluids 

The thermo-physical properties of nanofluid are func- 
tions of particle volumetric concentration and tempera- 
ture. The nanofluid density is given by [24,29] 

Pnf = 0Pnp + (1 - 0)Pwa 

where p np and p w is the nanoparticle density and 
water density, respectively, while cp is the particle volu- 
metric concentration. The nanofluid viscosity is esti- 
mated by [24] 

/x nf = £i exp (£ 2 </>)Mw, (5) 
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Figure 1 Schematic representation of (a) straight tube, (b) conical spiral tube, (c) in-plane spiral tube, and (d) helical spiral tube. 



where €\ and <£ 2 are constants (summarized in Table 
1), and [t w is the viscosity of base-fluid. 

The specific heat of nanofluid is assumed to be a weighted 
average of the base-fluid and the nanoparticles, e.g., 



Cp,nf - 



0PnpC P/ np + (1 - 0)PwCp, v 
Pnf 



(6) 



where c p>np and c p>w are the specific heats of nanopar- 
ticle and water, respectively. In this model, the thermal 
conductivity considers a combination of the static part 
of Maxwell's theory and the dynamic part taking the 
contribution of the Brownian motion of nanoparticles, 
defined as [24] 



fenf = 



fe np + 2fe w - 2(fe w - fe np ) (j) 
fe np + 2fe w + (fe w - fe np ) 0 



fe w + tifi(/)p w C piV 



KT 



-f{T,<p), (7) 



where d np is the nanoparticle diameter, ti is the Brow- 
nian motion constant, /c np and k w are thermal conductiv- 
ity of nanoparticle and water, respectively. Here, the 
effect of temperature and particle volumetric concentra- 
tion is taken into account in the Brownian motion from 
empirical data given by [24] 



0 = 0i(lOO0)* 



(8) 



/(T, 0) = (ci0 + c 2 )T/T 0 + (c 3 0 + c 4 ), (9) 
where /3 lt f} 2 > c h c 2, ^3 and C4, are constants (see Table 1). 



Thermo-physical properties of base-fluids 

The base-fluid considered in this article is water. 
Thermo-physical properties of water were obtained as 
polynomial functions of temperature [43]; the water 
density is defined by 



p w = -3.570 x KT 3 T 2 + 1.88T + 753.2, 
while the water viscosity is given by 

238.3 



/x w = 2.591 x KT 5 x IOT- 143.2 ; 



(10) 



(11) 



and the thermal conductivity of water is calculated 
from 

fe w = -8.354 x KT 6 T 2 + 6.53 x KT 3 T - 0.5981.(12) 
The specific heat of water is considered constant at 

c PiW = 4200. (13) 
Properties of nanoparticles are given in Table 1. 

Heat transfer performance 

The heat transfer performance of the cooling channel is 
discussed in terms of the figure of merit, FoM, which is 
defined as 



FoM 



W 



(14) 



pump 
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Table 1 Base case and operating parameters 



r al dlllclci 


Value 


Unit 


C p,np,AI203 


/ OJ 


j Ky ■ i\ 


^p,np, CuO 




j • Kg ■ i\ 


A 

"np, AI203 


J7 X I U 


m 


A 

"np, CuO 


1Q v 1 rf^ 
zy x i u 


m 


K np, AI203 


36 


\A/ v-W^ V~ 

vv ■ m ■ K 


L 

K np, CuO 


1 Q 
I O 


\A/ rv^~^ V~ 

vv ■ m ■ k 


ft 


j X I U 




ft 


I Jo I X I U 




Pnp, AI203 


ooUU 


kg ■ m 3 


Pnp, CuO 


CO I U 


kg ■ m 3 


m in 


9 X 1 0" 3 


kg ■ s 


Pout 


1 01 325 


D-i 

ra 


T 
'0 


1QQ 1 c; 
Zyo. I j 




T 
'in 


298.1 5 


I/ 


T 

'wall 


323.1 5 


1/ 


Cl 


Z.oZ 1 / X 1 U 




^2 


3 nil w 1 n _ 3 

o.y 1 / x 1 u 


_ 


c 3 


-o.uooy X I U 


_ 


C 4 


-3.91 123 x 10" 3 




Ci (Al 2 0 3 ) 


0.9830 




£2 (AI2O3) 


12.959 




£l (CuO) 


0.9197 




£ 2 (CuO) 


22.8539 




ft (Al 2 0 3 ) 


8.4407 
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where Wp ump is the pumping power required to drive 
the fluid flow through the channel. It is given by 

1 

W pump = mAp. (15) 

^pump 

Here, 7]p ump is the pump efficiency (assumed to be 
70%), W is the total heat transfer rate, and Ap is the 
pressure drop in the cooling channel. The total heat 
transfer rate is given as 

W = mc P/n f(T m/in - T m ,out)/ (16) 

where m is the mass flow rate and T m>in and r m>out are 
mixed mean temperature at the inlet and outlet, respec- 
tively. The mixed mean temperatures is calculated as 

Tm = ^vf TudAc ' (17) 

where A c is the cross section area of the channel and 
V is the mean velocity given by 



= }-J*dA c . (18) 



Boundary conditions 

The boundary conditions for the flow inside the channel 
are prescribed as follows 

♦ Inlet At the inlet, we prescribe inlet mass flow rate 
and inlet temperature. 

m = m m , T = T in . (19) 

♦ Outlet At the outlet, we specify the pressure and 
streamwise gradient of the temperature is set to 
zero; the outlet velocity is not known a priori but 
needs to be iterated from the neighboring computa- 
tional cells. 

p = pom, n • (fenfVT) = 0. (20) 



♦ Walls At walls, we set no slip condition for veloci- 
ties and constant wall temperature. 

u = 0, T = T wall . (21) 

In this article, a constant mass flow rate at a Reynolds 
number (Re = pUD h l { \a) of approximately 1000 is pre- 
scribed at the inlet for comparison purposes. 

Numerics 

The computational domains (see Figure 2) were created 
in AutoCAD 2010; the commercial pre-processor soft- 
ware GAMBIT 2.3.16 was used for meshing, labeling 
boundary conditions and determines the computational 
domain. Three different meshes, 1 x 10 5 , 2 x 10 5 , and 4 
xlO 5 , were tested and compared in terms of the local 
pressure, velocities, and temperature to ensure a mesh 
independent solution. It is found that mesh number of 
around 2 x 10 5 gives about 1% deviation compared to 
mesh size of 4 x 10 5 ; whereas the results from mesh 
number of 1 x 10 5 deviate by up to 8% compared to 
those from the finest one. Therefore, a mesh of around 
2 x 10 5 (20 x 20 x 500) elements was considered suffi- 
cient for the numerical investigation purposes; a fine 
structured mesh near the wall to resolve the boundary 
layer and an increasingly coarser mesh in the middle of 
the channel to reduce the computational cost. 

Equations 1-3 together with appropriate boundary con- 
ditions and constitutive relations comprising of five 
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dependent variables, u, v, w, p, and T, were solved using 
the finite volume solver Fluent 6.3.26. User-Defined func- 
tions (UDF) were written in C language to account for 
particle volumetric concentration and temperature-depen- 
dence of the ther mo -physical properties of the nanofluids. 

The equations were solved with the well-known Semi- 
Implicit Pressure-Linked Equation (SIMPLE) algorithm, 
first-order upwind discretization and Algebraic Multi-grid 
(AMG) method. As an indication of the computational 
cost, it is noted that on average, around 200-500 iterations 
and 500 MB of Random Access Memory (RAM) are 
needed for convergence criteria for all relative residuals of 
10" 6 , this takes 5-30 min on a workstation with a quad- 
core processor (1.83 GHz) and 8 GB of RAM. 

Results and discussion 

The numerical simulations were carried out for four dif- 
ferent tube geometries, four different nanofluid concen- 
trations, and two different nanofluid suspensions. The 
base-case conditions together with the physical para- 
meters are listed in Table 1, while the geometric details 
can be found in Table 2. 

Validation 

When developing and implementing mathematical model to 
predict the behavior of nanofluid heat transfer, one needs to 
pay special attention to validation of the model due to 



inherent complexity of coupled physical phenomena and 
interaction between base-fluid and nanoparticle. In this study, 
we aim to validate our model with an experimental nanofluid 
heat transfer by Anoop et al. [19], which has error of approxi- 
mately 4%. The heat transfer performance of nanofluid flows 
in circular tube with diameter 4.75 x 10" 3 m and length of 1.2 
m is approximated with 2D axisymmetric model, see Anoop 
et al. [19] for details of the experimental setup. 

The validation is initiated with heat transfer perfor- 
mance of water flowing at a constant Reynolds approxi- 
mately 1580; after which, the heat transfer performance of 
4 wt% of water- A1 2 0 3 nanofluid with nanoparticle size 45 
nm flows at Reynolds approximately 1588 is compared, as 
depicted in Figure 3. It is found that the model predictions 
agree well with the heat transfer performance from 

Table 2 Geometric parameters 
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experimental counterpart for both water and nanofluid. 
This implies that the model correctly accounts for the fun- 
damental physics associated with nanofluid heat transfer. 

Effect of geometry 
Base-fluid 

One of the key factors that determine the heat transfer per- 
formance is the cross-sectional tube geometry. This study 
examines four different square cross section tubes geome- 
tries: straight, conical spiral, in-plane spiral, and helical 
spiral with water as the base working fluid. Since the con- 
vective heat transfer inside the tube is directly linked to 
flow behavior, it is of interest to investigate the flow pat- 
terns inside the tubes. In our previous studies [1-3], albeit 
using air as working fluid, showed that the presence of cen- 
trifugal force due to curvature leads to significant radial 
pressure gradients in flow core region. In the proximity the 
inner and outer walls of the coils, however, the axial velo- 
city and the centrifugal force will approach zero. Hence, to 
balance the momentum transport, secondary flow should 
develop along the outer wall. This is indeed the case, as 
can be seen in Figure 4, where the secondary flow with 
higher velocities is generated in the outer wall region of 



coiled tubes (see Figure 4b,c,d). However, this is not the 
case for the straight tube (Figure 4a) as a fully developed 
flow exists inside the tube. It is noted that at this particular 
Reynolds number (approximately 1000), the secondary 
flows appear as one-pair for conical spiral and helical spiral 
tubes; whereas in the in-plane spiral tube, the secondary 
flows appeared as two-pairs. 

The presence of secondary flow with high velocities is 
expected to have direct impact on the heat transfer rate. 
This can be inferred from Figure 5 which presents tem- 
perature distribution over the cross sections of various 
tube designs. As can be seen from Figure 5, temperatures 
in coiled tubes are higher than in straight tube at the 
same axial distance which indicates that coiled tubes 
have higher heat transfer rate when compared to that of 
the straight tube due to the presence of secondary flows. 
It is also worth noting that the higher intensity of sec- 
ondary flow will tend to lead to higher heat transfer rate. 

Now looking at the mixed mean temperature and total 
heat transfer variation along the tube length (see dotted 
line in Figure 6), it is noted that coiled tubes have superior 
heat transfer performance when compared to that of the 
straight tube; the total heat transfer rate can be up to 




o 1 1 1 1 1 1 

0 50 100 150 200 250 

x/D 

Figure 3 Comparison of heat transfer coefficient between simulation (lines) and experimental data [19] (symbols) for water and 
nanofluid. 
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Figure 4 Velocity profiles of water flow in (a) straight duct; (b) conical spiral duct; (c) in-plane spiral duct; and (d) helical spiral duct at 
L = 50 cm. 



almost three times higher than that for the straight tube. 
In the near-inlet region, the heat transfer performance of 
in-plane spiral yields the best result among others, fol- 
lowed by conical spiral and helical spiral; whereas, in the 
near-outlet region, the helical coil performs the best fol- 
lowed by in-plane spiral and conical spiral This indicates 
that, for water as working fluid, in-plane spiral is more 
effective to be used in short tube applications, while the 
helical spiral is more effective for long tube applications in 
terms of amount of heat transferred. 
Nanofluids 

Four square cross section tube geometries were examined 
for flow of nanofluid suspensions of water- A1 2 0 3 with 
nanoparticle concentration of 1%. The results are depicted 
in Figure 6 where the mixed mean temperature and total 
heat transfer of base-fluid and nanofluids are shown. It is 



noted that adding 1% concentration of A1 2 0 3 in water 
improves the heat transfer performance. The total heat 
transfer for straight tube increases up to 50% as compared 
to that for water, whereas for coiled tubes, the heat transfer 
improves by about 50% in the near-inlet region and then 
decreases toward the outlet. Furthermore, among the coiled 
tube geometries, in-plane spiral gives the highest heat 
transfer improvement, followed by helical spiral and conical 
spiral tubes. This implies that in-plane spiral tube may have 
potential application to be used along with nanofluid due 
to its higher heat transfer performance. Therefore, the most 
of the following results refer to in-plane spiral coils. 

Effect of nanoparticle concentration 

The amount of nanoparticles suspended in the base- 
fluid plays a significant role in deter-mining heat 
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Figure 5 Temperature distribution of water flow in (a) straight duct; (b) conical spiral duct; (c) in-plane spiral duct; and (d) helical 
spiral duct at L = 50 cm. 



transfer performance. Intuitively, adding larger amount 
of nanoparticles in the base-fluid increases thermal 
conductivity of the nanofluid; however, care has to be 
taken as it also increases the friction factor and may 
reduce the stability of nanofluids due to agglomeration 
and sedimentation. To study the impact of these fac- 
tors, we investigated four different nanoparticle con- 
centrations: 0, 1, 2, and 3% of A1 2 0 3 in the base-fluid 
(water). Figure 7 displays the velocity profiles for the 
in-plane spiral tube for various nanoparticle concentra- 
tions. Interestingly, the velocity profiles are not 
strongly affected by the additional nanoparticle suspen- 
sion, especially at low concentrations. We note that at 
1 and 2% of A1 2 0 3 concentration, there is no signifi- 
cant difference on the secondary flow development 
inside the tube; whereas, at 3% A1 2 0 3 concentration, 



the effect of nanofluid suspension becomes stronger: 
the secondary flow appears in two-pairs as compared 
to that in one-pair at lower nanoparticle concentra- 
tions. A plausible explanation is the fact that nanofluid 
suspension does not significantly change viscosity of 
the fluid. Conversely, this is not the case for thermal 
conductivity of the nanofluid, as mirrored in Figure 8, 
where the addition of small amount of nanoparticle 
(1%) drastically changes the temperature profiles inside 
the tube. Furthermore, the temperature profiles for 
higher amount of nanoparticle concentration (2 and 
3%) also slightly change, but they are mainly affected 
by the hydrodynamics (secondary flows). 

Proceeding to the local mixed mean temperature and 
total heat transfer along the tube, as illustrated in 
Figure 9, it is clearly seen that additional small amounts 
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Figure 7 Velocity profiles of (a) water, (b) water with 1% Al 2 0 3/ (c) water with 2% Al 2 0 3/ and (d) water with 3% Al 2 0 3 flows inside an 
in-plane coiled tube at L = 50 cm. 
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Figure 8 Temperature distribution of (a) water, (b) water with 1% Al 2 0 3/ (c) water with 2% Al 2 0 3/ and (d) water with 3% Al 2 0 3 flows 
inside an in-plane coiled tube at L = 50 cm. 
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Figure 9 (a) Mixed mean temperature and (b) total heat transfer at various concentrations of Al 2 0 3 inside an in-plane coiled tube 
along the tube length. 
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of nanoparticles improves the heat transfer performance 
significantly, especially in the near-inlet area. How-ever, 
increase in nanoparticle concentration leads to a reduc- 
tion of total heat transfer along the tube by approxi- 
mately 5%. It is noteworthy that adding large amounts 
of nanoparticles in the suspension is not effective in 
enhancing heat transfer. Moreover, low nanoparticle 
concentration also has advantages of better stability of 
the suspension as it minimizes agglomeration and 
sedimentation. 

Effect of nanofluid type 

So far, the simulated nanofluid type chosen was water- 
A1 2 0 3 ; it is, therefore, of interest to see the heat transfer 
performance for a different nanofluid. In this study, we 
compare the performance of water-Al 2 0 3 and water- 
CuO nanofluids. Note that other types of nanofluid 
suspensions can be easily simulated within the frame- 
work of this model once their properties are known. 
Figure 10 shows temperature profiles for an in-plane 
spiral tube flowing through with water (Figure 10a), 1% 
of A1 2 0 3 nanofluid (Figure 10b) and 1% of CuO nano- 
fluid (Figure 10c). We note that the temperature profiles 
for both nanofluids (Figure 10b,c) are much higher than 
that of water (Figure 10a). Closer inspection reveals that 
a slightly larger area of higher temperature exists for the 
A1 2 0 3 suspension (Figure 10b) as compared to that for 
CuO suspension (Figure 10c). This is attributed to the 
stronger secondary flow observed in A1 2 0 3 nanofluid 
when compared to that of the CuO nanofluid (not 
shown here due to page limitation). 

The heat transfer performance of two different nano- 
fluid types is further evaluated in terms of the local 
mixed mean temperature and total heat transfer. As 
seen in Figure 11, the mixed mean temperature for the 
nanofluid is around 15% higher than that of water. 
There is no discernible difference between A1 2 0 3 and 
CuO suspensions in terms of the mixed mean tempera- 
ture. For total heat transfer, A1 2 0 3 gives somewhat 
higher heat transfer (approximately 5%) when compared 
to the CuO nanofluid. Therefore, it can be deduced that 
A1 2 0 3 nanofluid performs better heat transfer perfor- 
mance than that of CuO nanofluid, but not significantly. 
The stability and cost would decide the selection 
between these two nanofluids. 

Overall heat transfer performance 

A summary of heat transfer performance for all cases 
considered in this article is presented in Figure 12. Here 
several features are apparent; foremost among them is 
that the coiled tubes provide significantly higher heat 
transfer than that of straight tube, and addition of a 
small amount of nanoparticles in the base-fluid 
enhances heat transfer further (see Figure 12a). It is 
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Figure 10 Temperature distribution of (a) water, (b) water with 
1% Al 2 0 3/ and (c) water with 1% CuO flows inside an in-plane 
coiled tube at L = 50 cm. 
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Figure 12 (a) Total heat transfer, (b) pressure drop, and (c) Figure of Merit (FoM) of water and nanofluids in various tubes. 
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noted that the maximum heat transfer performance is 
achieved at 1% nanoparticle concentration, decreasing 
with higher amounts of nanoparticles. 

Aside from higher heat transfer performance, keeping 
pressure drop at a minimum is of interest for reducing 
the operating cost and saving energy. Figure 12b shows 
a summary of the pressure drop required for all cases 
studied. Note that the mass flow rate is kept constant in 
all cases; hence, it can be used directly to represent the 
pumping power required. The straight channel requires 
the lowest pressure drop among all cases; whereas the 
coiled tube designs require more than double the pres- 
sure drop of the straight channel. Among the coiled 
tubes, helical spiral tube needs the highest pressure 
drop, followed by in-plane spiral and conical spiral 
tubes. An interesting phenomenon is observed at a 
nanofluid concentration of 1% when the pressure drop 
for coiled tubes is slightly lower than that for water. 
This is due to the fact that at low particle concentra- 
tions, the particle volumetric concentration affects the 
nanofluid viscosity negligibly while the effect of tem- 
perature increases in the nanofluid thermo-physical 
properties. 

With respect to the heat transfer performance and 
pressure drop required in the system, the "Figure of 
Merit" concept is introduced as a measure of the heat 
transferred per unit pumping power (see Equation 14 
for details). Figure 12c presents the computed figures of 
merit for various tube geometries, nanofluid concentra- 
tions and nanofluid types. It is found that apart from 
the higher heat transfer rate, the coiled tubes have lower 
figures of merit than those of the straight tube. This can 
be explained by the higher pressure drops required in 
the coiled systems (see Figure 12b). Among all coiled 
tubes tested, the conical spiral tube gives the highest fig- 
ure of merit, followed by in-plane spiral and conical 
spiral tubes. Furthermore, for the straight tube, the addi- 
tion of nanoparticles improves the figure of merit signif- 
icantly, albeit it decreases with increasing concentration. 
For coiled tubes filled with nanofluids, on the other 
hand, the improvement of figure of merit is only shown 
at low particle concentration of 1% and then it drops 
lower than that of water when more nanoparticles are 
added. Clearly, these results suggest that one can add 
nanoparticle up to 1% volumetric concentration to 
water to enhance heat transfer performance in coiled 
tubes; higher nanoparticle concentrations are not 
recommended. 

Concluding remarks 

A computational study was conducted to investigate the 
laminar flow heat transfer performance of square cross 
section tubes, i.e., straight, conical spiral, in-plane spiral, 
and helical spiral, with water and two nanofluids. It is 



found that adding 1% nanoparticle volumetric concen- 
tration improves heat transfer performance and the fig- 
ure of merit for all tubes. However, higher amounts of 
nanoparticles is not recommended. In-plane spiral tubes 
give better performance than other coiled tubes for 
nanofluids. Furthermore, A1 2 0 3 nanofluid gives slightly 
better heat transfer performance than CuO nanofluid in 
coiled tubes. Future study will evaluate various modeling 
approaches for nanofluid heat transfer, e.g., single-phase, 
two-phase mixture, Euler-Euler, and Euler-Lagrange 
models, in coils with respect to the effect of secondary 
flow to the nanoparticle concentration. 
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